Uniform and staggered magnetizations induced by Dzyaloshinskii-Moriya interactions 
in isolated and coupled spin 1/2 dimers in a magnetic field 
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We investigate the interplay of Dzyaloshinskii-Moriya interactions and an external field in spin 
1/2 dimers. For isolated dimers and at low field, we derive simple expressions for the staggered 
and uniform magnetizations which show that the orientation of the uniform magnetization can 
deviate significantly from that of the external field. In fact, in the limit where the D vector of 
the Dzyaloshinskii-Moriya interaction is parallel to the external field, the uniform magnetization 
actually becomes perpendicular to the field. For larger fields, we show that the staggered magneti- 
zation of an isolated dimer has a maximum close to one-half the polarization, with a large maximal 
value of 0.35 g[is in the limit of very small Dzyaloshinskii-Moriya interaction. We investigate the 
effect of inter-dimer coupling in the context of ladders with Density Matrix Renormalization Group 
(DMRG) calculations and show that, as long as the values of the Dzyaloshinskii-Moriya and of the 
exchange interaction are compatible with respect to the development of a staggered magnetization, 
the simple picture that emerges for isolated dimers is also valid for weakly coupled dimers with 
minor modifications. The results are compared with torque measurements on Cu2(C5Hi2N2)2Ci4. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.40.Mg, 75.30.Kz 
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I. INTRODUCTION 



In Mott insulators, the Heisenberg interaction JS,;.Sj 
is in most cases the dominant source of coupling between 
local moments, and most theoretical investigations are 
based on modeling in which only this type of interac- 
tion is included. It has been known for a very long 
time, however, that other, less symmetric, interactions 
are present. For instance, unless there is an inversion 
center on a bond, spin-orbit coupling induces an anti- 
symmetric interaction of the form D ■ [Si x Sj), which is 
known as the Dzyaloshinskii-Moriya (DM) interaction.^ 
Since it breaks the fundamental SU(2) symmetry of the 
Heisenberg interactions, the DM interaction is at the 
origin of many deviations from pure Heisenberg behav- 
ior, such as canting 3 , or small g a p Sl | i S 4 SiL2i2ii£ It is also 
known to have a dramatic impact on the properties of 
antifcrromagnets in a magnetic field. Numerous experi- 
mental investigations of quantum antiferromagnets cur- 
rently in progress in large field facilities call for a de- 
tailed understanding of this problem jiiiSiiiiiiiSi Several 
issues have recently been the subject of rather inten- 
sive research. For instance, the impact on triplon Bosc- 
Einstein condensationi§iiLi& of DM interactions has been 
analyzed^ The interplay of frustration and DM interac- 
tions has also received significant attention. 21,22,23 The 
consequence of the breaking of SU(2) symmetry on the 
excitation spectrum is also well understood thanks to 



the work of several people including some of the present 
authors 4i*iSi2i2i2»i2i It is by now well established that a 
DM interaction can open a gap in otherwise gapless re- 
gions. The scaling of this gap with the magnitude of the 
DM interaction has been worked out for several cases&S 

Surprisingly, however, the other important conse- 
quence of the breaking of the SU(2) symmetry on the 
ground state properties of weakly coupled dimers, namely 
the development of a local magnetization, has not re- 
ceived much attention so far, although it is of immediate 
relevance to several compounds. It was shown in the case 
of SrCu2(BC>3)2 that a DM interaction can lead to the 
development of a measurable (and in fact quite large) 
staggered magnetization^ but a simple picture of how 
the magnitude and the orientation of the DM interac- 
tion with respect to the magnetic field influences these 
properties has not yet emerged. Besides, the fact that a 
DM interaction can lead to the development of a trans- 
verse uniform magnetization and its impact on torque 
measurements of the magnetization have not been in- 
vestigated in detail. All these questions are central to 
the understanding of several systems of current inter- 
est. In particular, recent NMR results by Clemancey 
et alisl have revealed the presence of a staggered magne- 
tization in the dimer compound Cu2(C5Hi2N2)2Ci4 [ab- 
breviation: Cu(Hp)Cl], and the interpretation of these 
results requires a precise investigation of the effect of DM 
interactions on weakly coupled dimer systems. 

In this paper, we present a systematic analysis of the 
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consequences of an intra-dimer DM interaction on the 
development of local magnetization in systems of weakly 
coupled spin 1/2 dimers. We first look at the case of an 
isolated dimer, and derive simple expressions in the lim- 
its of weak and strong magnetic field which we believe 
are very useful to get a simple picture of subtle effects 
such as the effect of the relative orientation of the mag- 
netic field and the D vector of the DM interaction on 
the uniform magnetization. We then turn to the case of 
coupled dimers and concentrate on a simple ladder ge- 
ometry. This choice is motivated partly by the potential 
relevance of this geometry to actual compounds such as 
Cu(Hp)Cl, and by the possibility to obtain very accurate 
results using the Density Matrix Renormalization Group 
method (DMRG)2£i2£i21 in this quasi-one dimensional ge- 
ometry. Finally, we report new torque measurements on 
Cu(Hp)Cl and discuss them in the light of these results. 



II. ISOLATED DIMER 

The problem of an isolated dimer in a magnetic field 
in the presence of a DM interaction is defined by the 
Hamiltonian 

H = JSi ■ S 2 + D • (Si x S 2 ) - gfi B H(S z i + S z 2 ). (f) 

The z axis has been chosen to be that of the mag- 
netic field, and the yz plane as the plane defined by the 
magnetic field and the D vector (see Fig. P). In actual 
systems, the direction of the D vector is fixed by the 
microscopic arrangement of atoms and orbitals, and it 
is the orientation of the magnetic field that can be var- 
ied with respect to the crystal, but the convention of 
having the magnetic field along the z axis makes the dis- 
cussion somewhat simpler. The D vector is written as 
D = (0, D sin 6>, D cos 0). 
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FIG. 1: Pictorial representation of the model of Eq. Q of a 
dimer with DM interaction in a magnetic field. 

The isolated dimer problem is, of course, very simple. 
The Hilbert space is of dimension 4, and it will prove 
convenient to work in the basis 

\s) = ^(ITD-UT)), 



l*o> = ^=(IU) 
\t-i) = I u>. 



IT)), 



(2) 



Unfortunately, the explicit expressions that can be de- 
rived for the eigenenergies and eigenvectors are cumber- 
some and not particularly useful. Therefore, we will dis- 
cuss the results from the point of view of symmetry, will 
derive useful expressions for small D in weak field and 
close to saturation, and will present plots of some repre- 
sentative exact results. 



A. Symmetry analysis 

In a magnetic field, the SU(2) symmetry of the Heisen- 
berg model is reduced to a U(l) symmetry corresponding 
to a rotation around the field direction. 

As soon as a DM interaction with a D vector not par- 
allel to the magnetic field is introduced, the rotational 
symmetry in spin space is completely lost. The only sym- 
metry that remains is the mirror symmetry with respect 
to the yz plane (the plane containing the magnetic field 
and the D vector), which exchanges sites 1 and 2 and si- 
multaneously changes the sign of the x-component of the 
spin operators (the component perpendicular to the mir- 
ror plane). As a consequence, the expectation values of 
local spin operators in any eigenstate of the Hamiltonian 
satisfy the relations: 



(Si) 
(Sf) 



(SI) 

(S z 2 ) ■ 



(3) 



These relations imply that the staggered magnetization 
per site, defined as m s = ((Si — Sa))/2, is perpendicular 
to the plane defined by the magnetic field and the D 
vector, while the uniform magnetization per site defined 
by m u — ((Si + S2))/2 must lie in that plane. 

If the D vector is parallel to the field, the U(l) ro- 
tational symmetry is still present, which can be easily 
checked since S l ° l = S* + S% commutes with SfS% — 
5fS*f . The states and \ti) are still eigenstates with 
energies J/4 ± g^sH, but \s) and \to) get coupled. The 
staggered magnetization is identically zero, while the uni- 
form magnetization jumps abruptly from to 2g/iBZ (z 
is the direction of the applied magnetic field) at a critical 
field H c larger than its D — value J/gfiB- 



B. Low- field limit 

In the limit D/J <C I and below the saturation field 
H c = J/gfig, the ground-state wave function, up to sec- 
ond order in D/J, reads: 



TT), 



J>sin6> 



1 ' 4j0 |s> 2 v 2fJ nnuih 



h) 
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D sin 9 



2a/2( J + 9HbH) 



!*-!>• (4) 



In the low- field limit, first-order perturbation theory 
in H can be used to derive simple expressions for the 
expectation value of the various spin operators: 



(Sf) = -<Sf) = 

(sf) = (s%) = 

(Sf) = (SI) = 



g^sHD sin 9 
2J2 ' 
g/j, bHD 2 cos 9 sin t 
4J3 

g^sHD 2 sin 2 
4J3 ' 



(5) 



These expressions lead to compact and suggestive expres- 
sions for the uniform and staggered magnetizations: 



m. 



4J 3 
ffMg 
2J 2 



(D x H) x D, 
(D x H) . 



(6) 



As required by symmetry, the staggered magnetization 
is perpendicular to both the field and the D vector. As 
far as the uniform magnetization is concerned, symmetry 
only requires that it lies in the plane of the magnetic field 
and of the D vector, but in the low field limit, Eq. © 
shows that it is perpendicular to the D vector. So the 
uniform magnetization is in general not parallel to the 
magnetic field, as it would be in a system with SU(2) 
symmetry, and it can in fact deviate strongly: In the 
limit where the D vector becomes parallel to the field, 
the uniform magnetization becomes perpendicular to the 
magnetic field, a rather anomalous behavior that should 
have important consequences for torque measurements of 
the magnetization. 

Another remarkable feature of these results is that the 
staggered magnetization is first order in D, while the uni- 
form magnetization is second order. Thus, at low field 
the response is dominated by the staggered magnetiza- 
tion, as already observed in SrCu(B03)2- 

Finally, let us emphasize that, as implied by Eq. 
the uniform and staggered magnetizations have universal 
expressions in terms of the magnetic field and of the D 
vector, which are valid regardless of their orientation with 
respect to the lattice. 
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FIG. 2: Field dependence of the uniform and staggered mag- 
netizations per site m^, ra£ and mj for D/J = 0.04 and 
9 — 7r/4. The dashed lines are the analytical results derived 
in the limit D/J <C 1. Note the difference of scale for positive 
and negative magnetizations. 



large in the limit where D goes to zero. Note, however, 
that the staggered magnetization only takes significant 
values close to H = J/gfi b, in an interval of width of the 
order of D sin 9 which shrinks to zero in the limit where 
D goes to zero, consistent with a vanishing staggered 
magnetization when D = 0. 

When D -C J, the uniform magnetization at this field 
is equal to gfiB > which corresponds to half the polariza- 
tion value. When the angle between D and H is not 
7r/2, a small uniform component develops along y due 
to the coupling of \s) with \to). This transverse (with 
respect to the field) uniform magnetization is given by 

m v = -(v / 2/4)cos9{D/J)gfi B ^ -0.35 cos9(D/J)gn B - 
In contrast to the small field result, it is now linear in D, 
but remains much smaller than the staggered magnetiza- 
tion, which is of order one. 



C. Critical field 

At the critical field H c — J/g^s, one has to turn to 
degenerate perturbation theory since, for D — 0, \s) and 
\ti) are degenerate. When the D vector is not parallel 
to the field, these states get coupled by an off-diagonal 
term D sin 9. The ground state wave function is then 
simply given by (f>o — (\s) — \ti))/y/2, and the staggered 
magnetization per site is equal to (y/2/4)gfj, B ~ 0.35 g\is- 
Interestingly enough, this maximal value is independent 
of the angle 9 and does not vanish but remains quite 



D. Exact results 

To get an idea of the accuracy of the expressions ob- 
tained at low field and close to the saturation field, we 
have plotted in Fig. [21 the exact value of mf, m\ and 
m z u for a representative case (D/J = 0.04 and 9 = 7r/4). 
The small field expression is quantitatively accurate up 
to H ~ 0.25J/gfiB, and the width of the peak of the 
staggered magnetization and the maximal value of 
are indeed of order D. 
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III. COUPLED DIMERS (LADDER) 
A. The model 

In this section, our goal is to check to which extent the 
properties of a system of weakly coupled dimers resemble 
those of isolated dimers. In particular, the transition 
between zero magnetization and polarization takes place 
through an extended region of magnetic field of the order 
of the inter-dimer coupling, and we would like to know 
how the system behaves within and outside this region. 
We will attack this problem numerically, and in order to 
perform simulations on large systems, we have chosen to 
work in a ladder geometry and to use the DMRG. The 
model is defined by the Hamiltonian 

H = J^S M .S ii2 + ]r(-ir D-(S M xS M ) 

i i 

+ J\\ / J(S^i ■ Si+1,1 + 8i t 2 ■ St+1,2) 

i 

- g» B Hj2(S? A + Sl 2 ) (7) 

i 

As for the isolated dimer, the D vector is assumed to lie 
in the yz plane, i.e., D = (0, DsmO, Dcos8). Our choice 
of an alternating D vector from one rung to the other 
(see Fig. |3J) is motivated by symmetry considerations. 
Indeed, in a canonical ladder, the middle of each rung is 




FIG. 3: Ladder with staggered DM interaction. 

an inversion center, and the DM interaction vanishes by 
symmetry. A simple way to allow for the DM interaction 
to become finite without modifying the symmetry of the 
exchange couplings is to assume that some buckling is 
present along the ladder, as sketched in Fig. QJ In that 
case, the only mirror plane that contains a bond is the 
xz plane, and a DM interaction with a D vector paral- 
lel to y is allowed by symmetry. But, in this geometry, 
the presence of a C 2 axis (see Fig. 0J) implies that the D 
vector alternates from one rung to the other. The buck- 
ling realized in (^(CsH^^^CLt is slightly more subtle 



(successive rungs are connected by an inversion symme- 
try in the middle of a plaquette), but this symmetry also 
implies alternating D vectors. Note, however, that other 
ways of breaking the inversion symmetry of the rungs can 
lead to other arrangements of D vectors. 




FIG. 4: Structure of a buckled ladder. In such a ladder, 
a staggered DM interaction in the ^/-direction is allowed by 
symmetry. 

Another motivation to work with alternating D vectors 
is to keep the perturbation caused by the inter-dimer cou- 
pling as small as possible. In that respect, this choice is 
natural. Indeed, as we have seen in the previous section, 
the presence of a D vector on a rung induces a staggered 
magnetization. If the D vectors of neighboring rungs i 
and i + 1 are equal, the moments (S^i) and (Si+^i) will 
also be equal, which is in conflict with antifcrromagnetic 
inter- rung exchange interactions. If, on the contrary, the 
D vectors are opposite on neighboring rungs, the local 
moments will adopt configurations that are compatible 
with the exchange. 

B. Symmetry analysis 

With this choice of staggered D vectors, the model 
possesses the following symmetries: i) a yz mirror plane 
going through the middle of the rungs; ii) an inversion 
center in the center of the plaquette formed by two con- 
secutive rungs; iii) even translation symmetries along y. 
As long as these symmetries are not broken, the follow- 
ing relations between the expectation values of local spin 
operators on two neighboring rungs are expected to be 
satisfied: 

= (SI2) = (<S?+i,i) = (^+1,2) ■ (8) 

We thus define the staggered and uniform magnetizations 
per site as 

m s = (l/iV)^(-ir((S. 1 , 1 >-(S l , 2 » 

i 

m u = (l/iV)^((S i)1 ) + (S i , 2 )) , (9) 
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FIG. 5: Examples of the variation of the uniform magnetiza- 
tion along the z axis with the field. Inset: Plot of as a 
function of (D sin#/ J) 2//5 slightly below H c i, which confirms 
the scaling predicted in Ref. 
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FIG. 6: Lower panel: y component of the uniform magne- 
tization as a function of the field for 9 = 7r/6 and various 
values of D/J. Significant values appear between H c \ and 
H C 2, and far outside this interval as soon as D/J is not too 
small. Upper panel: Angle <f> u between the magnetic field and 
the uniform magnetization m u as a function of the field for 
9 — 7r/6 and several values of D/J. Note that 9 — cj> u goes to 
7r/2 in the low-field limit, in agreement with the prediction 
for an isolated dimer (Eq. EJ. 



where N is the total number of sites, and with the con- 
vention that the angle 8 is positive for i even. As in the 
isolated dimer case, the staggered magnetization m s is 
along the x axis, while the uniform magnetization m u 
lies in the yz plane. 



C. Uniform and staggered magnetizations 

Let us now turn to the discussion of the numerical re- 
sults we have obtained for the model of Eq. (J7|. We 
are interested in the regime D < Jj| < J. For D = 0, 
the model is a simple ladder in a field, and the proper- 
ties are well understood. There is, of course, no stag- 
gered magnetization because of the U(l) symmetry, and 
the uniform magnetization is parallel to the field for the 
same reason. It vanishes below a critical field H c i, takes 
off with a square-root singularity, and reaches saturation 
with another square-root singularity at a second critical 
field H C 2- The difference H C 2 — H c \ scales with Jii . Since, 
apart from this scaling, the properties depend very little 
on Jii , we quote results for a single value of Jii , and hav- 
ing in mind the compound Cu(Hp)Cl^i we have chosen 
Jii/ J = 0.2. For that ratio, the critical fields in the ab- 
sence of DM interactions are given by <?/^B#ci = 0.82 J 
and gfisH C 2 = 1.40 J. 

For the model with DM interaction, we have performed 
Exact Diagonalization (ED) 27 up to 20 sites (10 rungs), 
and DMRG calculations on ladders with up to 80 rungs. 
The results evolve smoothly with the size, and we only 
quote DMRG results obtained for 80-rung clusters (finite- 
size effects for the gap are discussed in the next section) . 
Note that in those calculations, S z is not a good quan- 
tum number. This is well known to reduce greatly the 
maximal size available to exact diagonalizations, but this 
also has an impact on the number of states we were able 
to keep during the DMRG runs. Here, we diagonalize 
(by the Davidson method) a matrix of size 4m 2 at each 
DMRG step. In a standard DMRG run where S z is a 
good quantum number, the matrix of the effective Hamil- 
tonian in the variational basis is block-diagonal, which 
can speed up the diagonalization by a factor of 10 or 
more. The memory needed is also larger at fixed m than 
for the standard DMRG. For those reasons, most of the 
calculations were done with up to m — 600 states kept 
during 5 sweeps, and only up to N = 80 sites. The dis- 
carded weight was of the order of 10~ 10 when we targeted 
two states to extract the gap, and of the order of 10~ 12 
or less when we targeted a single state to extract cor- 
relations. We also performed a few runs with m up to 
800 in order to confirm that the numerical data were well 
converged. 

The z component of the magnetization is displayed in 
Fig. El for several values of D and 8. It is reminiscent of 
that for D = 0; however, when D ^ 0, the magnetization 
develops as soon as the magnetic field is switched on, 
only reaching saturation asymptotically in the limit of 
infinite field. The square root singularities are removed. 
It was shown in Ref. that, at H c i, the magnetization 
should depend on the magnitude of the D vector as oc 
(Dsin#) 2 / 5 , in agreement with the present results (see 
the inset of Fig. |SJ . 

When 8 ^ tt/2 (i.e., D z ^ 0), a uniform magnetization 
along the y axis also develops, as in the isolated dimer 
case. Fig.[S]shows the magnetization along the y axis and 
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FIG. 7: Staggered magnetization as a function of the mag- 
netic field for several values of D/J and 9. Large values are 
achieved between H c \ and H C 2, and far outside this interval 
as soon as D/J is not too small. The value between H c i and 
H C 2 depends relatively weakly on D/J and 9, and is of the 
same order as the maximal value in the case of an isolated 
dimer (0.35 ff/Lts). In contrast, the value outside this interval 
depends very strongly on the magnitude of D sin 9. 



the angle 4> u between the uniform magnetization and the 
z axis as a function of the magnetic field for 9 — tt/6. At 
low field, the uniform magnetization is orthogonal to the 
DM vector, again as for an isolated dimer. The magneti- 
zation along y is maximal between the two critical fields. 
Its value in that range is clearly much smaller than the 
component along the field (<p u becomes very small near 
Hci), but this extra contribution to the uniform magne- 
tization will produce a torque that should be detectable 
experimentally given the very high sensitivity of torque 
measurements. 

The staggered magnetization along x exhibits a kind of 
plateau in the intermediate phase between H Cl and H C2 
(Fig. |7|) . Its magnitude inside the plateau is of the order 
of the maximal value of the isolated dimer (0.35 g^s), 
and it depends relatively weakly on D. In contrast, the 
extent of the tails outside this plateau region increases 
rapidly with D. Remarkably, the magnetization per spin 
along x is larger than along z up to H c \ and even slightly 
above. Note that the staggered magnetization depends 
essentially on the value of D y and is very weakly affected 
by the value of D z . 



D. Gap 

The effect of a SU(2) breaking interaction on a lad- 
der has been studied in Ref. jig. It strongly depends 
on the nature of the plateau phase. For the transition 
from the zero or full polarization to the gapless phase, 
the effective field theory is expected to be the same as 
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FIG. 8: Upper panel: Field dependence of the excitation 
gap A for J\\/J = 0.2 and several values of D/J: D/J = 
O.Ol(black), D/J = 0.02(red) D/J = 0.04(blue) D/J = 
0.08(green) and N = 80 (DMRG). Lower left panel: Scal- 
ing of the excitation gap as a function of 1/N for D/J=0.02. 
Lower right panel: Scaling of the gap as a function of D/ J for 
N=80 (see text). 



for the spin chain close to saturation, and the gaps at 
H c i and H C 2 should open as (D sin #) 4//5 , as shown in 
Ref. p| . This prediction clearly agrees with the results 
for 9 — ir/2 shown in Fig. [S] (lower right panel). Size 
effects are already very small for N = 80 sites, as can be 
seen in Fig. [S] (lower left panel) . Between the two critical 
fields, the gap is expected to remain finite. (The closing 
of the gap in Ref.hfi was caused by a breaking of the 
Z2 symmetry which does not occur here as there is no 
to = 1/2 plateau when D/J = 0). The effect of the z 
component of D is expected to be very small. This is 
also confirmed by our DMRG results (not shown). 



E. Torque measurement 

If the uniform magnetization is not parallel to the field, 
it induces a torque r on the system. Usually, such a 
torque is only present if the field is not along a high 
symmetry direction of the g-tensor. This is the basis of 
torque measurements of the magnetization. However, as 
shown above, a DM interaction can also induce a compo- 
nent of the magnetization perpendicular to the magnetic 
field, which should show up in torque experiments as 
an additional contribution. Interestingly enough, torque 
measurements on Cu(Hp)Cl indeed reveal the presence 
of such a contribution. Experiments were carried in a 
resistive magnet and r was measured up to 23 T at 410 
mK. The orientation of the crystal was adjusted so that 
t = at the highest values of H, as shown in the in- 
set of Fig. which fully cancels the contribution due to 
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the anisotropy of the g tensor. This orientation indeed 
corresponds to ffii [100]. In spite of this, a large addi- 
tional contribution shows up between the critical fields, 
and extends well outside the intermediate region. For 
comparison, the calculated component of the uniform 
magnetization perpendicular to the field of a ladder with 
J,, I J = 0.2, D y /J = 0.05 and D z /J = 0.086 (9 = tt/6) is 
depicted on the same plot, with scales adjusted to get the 
same value at H c \. The values of D y and J\\ are those 
used in Ref. to fit the staggered magnetization, while 
the results depend very little on D z up to an overall scale 
factor. The two curves are in good qualitative agreement, 
especially considering the fact that the only adjustable 
parameter is the overall scale factor. In order to go be- 
yond this qualitative agreement, it would be necessary to 
consider several additional effects. First of all, inelastic 
neutron scattering data have challenged the description 
of this system as a simple ladder ^ and further couplings 
(still note definitely identified) should presumably be in- 
cluded. In addition, there is a transition into a 3D or- 
dered phaseii between H cl and H c2 , and although the 
precise nature of the ordering is still unknown, it is very 
likely it will affect the uniform magnetization. Clearly, 
at the present stage, too little is known about these ad- 
ditional effects to be able to take them into account, and 
this is left for future investigation. 




FIG. 9: (Color online) Transverse uniform magnetization m y 
for Dy/J = 0.05, D z /J = 0.086, J\\/J = 0.2 (blue circles) 
and torque divided by field (experimental curve obtained on 
Cu(Hp)Cl) as a function of the field. Inset: torque measure- 
ment raw data (black) and torque divided by field (red). 



IV. CONCLUSIONS 



If spin 1/2 dimers are coupled in such a way that there 
is no inversion center at the middle of the bond, very sig- 
nificant modifications of the physics in a magnetic field 
have to be expected. Indeed, unless it is forbidden by 
symmetry, a DM interaction will always be present, and 
the analysis reported in this paper shows that even a tiny 
DM interaction can modify some aspects of the physics 
rather dramatically. This is especially true for the stag- 
gered magnetization, which immediately acquires large 
values in the intermediate phase where the system gets 
polarized, and which can take on significant values out- 
side this phase for physically relevant values of the DM 
interaction. This is also true for the uniform magnetiza- 
tion as soon as the D vector of the DM interaction and 
the field are neither parallel nor perpendicular. In that 
case, a component of the uniform magnetization perpen- 
dicular to the magnetic field appears, which can induce a 
measurable torque on the sample. This has been proven 
for an isolated dimer and for a ladder with staggered 
DM interactions, but these conclusions are expected to 
hold true for all coupled-dimer systems as long as the 
D vectors are arranged in such a way that there is no 
competition with Heisenberg exchange as far as the de- 
velopment of a staggered magnetization is concerned. It 
is our hope that these results will help understand some 
of the strange properties observed in coupled-dimer sys- 
tems. 
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